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The spectrum of twelve-flavor QCD has been studied in details by the LatKMI collaboration. In 
this proceeding we present our updated results for the spectrum obtained with the HISQ action at 
two lattice spacings, several volumes and fermion masses. In particular, we emphasize the exis¬ 
tence of a flavor-singlet scalar state parametrically light with respect to the rest of the spectrum, 
first reported in our paper [1], This feature is expected to be present for theories in the con¬ 
formal window, but the lattice calculation of such a state is difficult and requires noise-reduction 
techniques together with large statistics, in order to evaluate disconnected diagrams. Being able to 
provide a robust observed connection between a light flavor-singlet scalar and (near-)conformality 
is an important step towards observing a light composite Higgs boson in walking technicolor the¬ 
ories on the lattice [2], We also show updated results for the mass anomalous dimension y„, 
obtained from various spectral quantities, including the string tension, under the assumption that 
the theory is inside the conformal window. 
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1. Introduction 

Four years ago the LatKMI collaboration started a program of lattice simulations and related 
theoretical work to systematically study SU(3) gauge theories with increasing number Nf of mass¬ 
less flavors in the fundamental representation. Lattice simulations of such theories were aimed 
at the discovery of a near-conformal theory with a mass anomalous dimension y,„ ~ 1, a possible 
candidate for Walking Technicolor [3]. 

For large number of flavors, the fi function of SU(3) gauge theories displays an infrared con¬ 
formal fixed point in its perturbative expansion [4]. This happens for a critical number of fla¬ 
vors Nf > N" f , where Nj is close to Nf = 8 for 2-loop perturbation theory or Nf = 12 for ladder 
Schwinger-Dyson analysis. The perturbative treatment is not completely justified at such value of 
Nf because the IRFP becomes strongly coupled. Non-perturbative lattice simulations are a suited 
tool for the study of this fixed point. 

The SU(3) 12-flavor theory has been explored on the lattice by many different groups and 
using several different methods. LatKMI investigated signals of a IRFP in this theory by looking 
first at the spectrum of mesonic bound states and their decay constants, specifically the lightest 
pseudoscalar and vector state [5]. Subsequently, following the results of Ref. [6], LatKMI inves¬ 
tigated the gluonic spectrum, focusing on the scalar glueball state [7] and the string tension [8]. 
The encouraging results of Ref. [7] indicated the presence of a light flavor-singlet scalar state in 
the spectrum, lighter than the pseudoscalar one, and this state was carefully studied in another 
publication using fermionic operators as well as gluonic ones [1], 

In this proceeding we report updated results to all of the aforementioned papers by LatKMI. 
The results presented in this publication are to be considered preliminary. The lattice simulations 
performed include larger volumes, smaller fermion masses and increased statistics with respect 
to published results. We use the Highly Improved Staggered Quark (HISQ) action and the tree- 
level Symanzik gauge action at two bare gauge coupling constants (5 = 6 /g 2 = 3.7 and 4.0 and 
we simulate four physical volumes with L = 18, 24, 30, 36 with aspect ratio T/L = 4/3. Different 
bare quark masses are used, going from nif = 0.2 on the smallest L = 18 volume, to m/ = 0.03 on 
the largest L = 36 (m/ = 0.04 for /3 = 4.0). As a reference, the simulated pseudoscalar (would-be 
pion) masses are 0.24 < M n < 0.91. Some of the features of these simulations, are the good flavor 
symmetry realization of the HISQ action and the large number of configurations, the latter needed 
for a precise determination of gluonic correlators and fermionic disconnected diagrams [1]. 

2. Spectrum and finite-size scaling 

We first focus on the lightest flavor non-singlet mesonic resonances of the lattice theory. The 
pseudoscalar mass ( M n ), the pseudoscalar decay constant (F K ) and the vector mass (M p ) are ex¬ 
tracted from lattice Euclidean correlators of the corresponding staggered fermions operators. The 
standard techniques used are the same described in Ref. [5]. The goal is to look at this hadronic 
spectrum as a function of the bare fermion mass m/ to identify signals of a conformal fixed point. If 
the theory is inside the conformal window the leading mass dependence of the spectrum is dictated 
by hyperscaling [9] relations: 

t i 

M H ^m l f +ym , F h oc m j +ym . (2.1) 
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Therefore dimensionless ratios of masses and decay constants should approach a constant value 
when the theory reaches the chiral limit where the IRFP is present. Moreover, we should remind 
that the relations in Eq. (2.1) are valid in the infinite volume limit and when the mass is the only 
source of soft breaking of scale invariance. 

In Fig. 1 we plot the dependence of two ratios, F n /M K and M p /M K , on M n as a proxy for the 
fermion mass my. For decreasing M n the theory is approaching the chiral limit and the dimension¬ 
less ratios first increase before reaching a plateaux value. These plateaux are only visible if larger 
volumes are considered at small pseudoscalar masses and this is a clear signal of how finite-volume 
effects can influence the system. The results presented in Fig. 1 improve the similar ones in Fig.4 
of Ref. [5] by adding data points on larger volumes and smaller pseudoscalar masses. 
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Figure 1: (Left) The ratio of the pseudoscalar decay constant and pseudoscalar mass at two different values 
of the coupling constant. (Right) The ratio of the vector mass and the pseudoscalar mass at two different 
values of the coupling constant. In both cases the ratios approach an approximately constant value towards 
the chiral limit of the theory, signaling a common scaling function. 


One can turn finite-volume effects, which distort hyperscaling relations, into a powerful tool 
to look for infrared conformality, in the same way one uses finite size scaling techniques to study 
critical points of finite systems in Statistical Mechanics. At finite fermion mass and finite lattice 
size, the presence of a conformal IRFP can be studied using universal scaling relations for hadronic 
masses and decay constants: 

$=LM h = fit (jc) , = LF n = f F (x) , (2.2) 

where the subscript H refers to a hadron (7t or p here) and the scaling variable is 

x = Lmf^ . (2.3) 

In Ref. [5] the value of y m is estimated by quantifying the alignment of the data points for 
i %{L,mf) to a universal function. The scaling functions fu (x) and f F (x) have to reach the asymp¬ 
totic infinite-volume form fx{x) ~ x required by Eq. (2.1). Moreover, the value of y m is character¬ 
istic of the critical point and must be independent of the observable used in the scaling analysis. 
The analysis presented in Ref. [5] showed that a common value of y m could be extracted by looking 
at both M k and M p . However, the data for F K did hint at a y m value that was not universal. This 
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would weaken any possible claim about the existence of a conformal IRFP in the chiral limit of the 
theory. It is insightful hence to notice that including new data points on larger volumes drastically 
changes the situation. This can be seen by comparing the two panels of Fig. 2, corresponding to 
old published data [5] and newly available data presented for the first time in this proceeding. From 
the plot one can read of a value 0.4 < y m < 0.5 for the anomalous mass dimension at the IRFP. 
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Figure 2: (Left) Optimal value of y m with statistical and systematic errors from Ref. [5]. (Right) Updated 
values of y m for three different observables and two coupling constants using larger volumes at small fermion 
masses. The y m obtained from the finite-size scaling analysis of F n is now compatible with the one coming 
from other channels. 

A different approach that includes corrections to the finite size hyperscaling relations can help 
in taking into account volume effects and non-universal corrections. For example, one could try 
and fit the data <§ (L,m/) using different functions: 


% = co + ciLtnJ (1+7 "' ) 

(2.4) 

= Co + c x Ltnf^ + C2Lnif 

(2.5) 

= (co + ciLffjy^ 1+7 "^)(l+c 2 m“) . 

(2.6) 


In Eq. (2.5) the value of a is fixed to a = (3 — 2y,„)/(l + y n ), inspired by solution of the Schwinger- 
Dyson equation [10], or a = 2, regarded as a lattice discretization artefact. In Eq. (2.6), the addi¬ 
tional term °c mj is inspired by the work in Ref. [11], where corrections due to the (near-)marginal 
gauge coupling operator have been shown to be relevant in certain parameter regions of the simu¬ 
lations. By using fit forms above on large volume data ( LM K >8.5 and LF n > 2) the fitted value y m 
turns out to be consistent across different observables and coupling constants, as well as consistent 
with the previous analysis shown in Fig. 2 (Right). Including corrections due to near-marginal op¬ 
erators near the fixed point (ft) f 0) gives compatible results and it also improves the quality of the 
fits in the heavy mass region. 


3. Flavor-singlet scalar 

When a conformal theory is perturbed by a mass deformation as in our lattice simulations at 
fixed nif, a flavor-singlet scalar excitation in the spectrum may be parametrically light with respect 
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to the other resonances. This non-perturbative physical phenomenon was observed first on the lat¬ 
tice using gluonic interpolating operators in the SU(2) theory with two adjoint flavors [6], and then 
observed also in SU(3) with twelve fundamental flavors [7, 1] using, for the first time, fermionic 
correlators. We employ large computational resources to measure disconnected fermionic contri¬ 
butions that are necessary for the correct analysis of the flavor-singlet scalar correlator. 

In order to study the flavor-singlet scalar correlator and its ground state, which we call a in 
analogy to QCD, we employ the local fermionic bilinear operator 

3 

^ct(0 = , (3-1) 

i— 1 x 

where the index i runs through different staggered fermion species. Then we calculate its corre¬ 
lator, which is constructed by both the connected C(t) and vacuum-subtracted disconnected D(t) 
correlators, {@a{t)@a{ 0)) = 3 D(t)—C{t), where the factor in front of D(t) comes from the number 
of species. Following the analysis performed in Ref. [1] we extract an effective mass for the state 
in the scalar channel and we study its dependence on the bare fermion mass in order to check if 
hyperscaling takes place. Due to the higher computational cost we focus on a fixed lattice spacing 
jS = 4.0. 

The left panel of Fig. 3 shows that the lightest scalar state is lighter than the pseudoscalar one 
in the whole mass region explored. The fermion mass dependence is well described by Eq. (2.1) 
when y„, is fixed to the value extracted from the study of M n . The right panel of Fig. 3 instead shows 
the comparison with the scalar spectrum obtained from gluonic interpolation operators (following 
the methodology of Ref. [12]). A scalar glueball seems to exist in the spectrum, but it is at the same 
scale of the vector meson. 
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Figure 3: (Left) Spectrum of the flavor-singlet scalar obtained on different volumes from fermionic interpo¬ 
lating operators. The pseudoscalar spectrum is shown for reference. Only the filled symbols are included in 
the hyperscaling fit. (Right) The spectmm in the scalar channel obtained from gluonic interpolating opera¬ 
tors. M n , M p and M a are all shown for reference. 


4. String tension 

We also study purely gluonic quantities to extend the hyperscaling analysis to more observ- 
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ables, which are affected by different systematic effects compared to the fermionic ones [8]. Among 
the gluonic observables explored in the SU(3) Nf = 12 theory is the string tension fG. Using two 
different type of correlators, we obtain compatible results for the string tension at /3 = 4.0 and for 
several fermion masses and volumes. The first method to obtain yfa is from smeared Polyakov 
loop correlators coupling to torelon excitations whose mass grows with the length of the loop 


nZ(or (T) — 6 L 


n 
3 L 


7t 2 1 

18 z 7 a ' 


(4.1) 


The second method uses APE smeared Wilson loop operators to construct Creutz ratios which 
describe the static quark-antiquark potential 


V(r) =v 0 - y + (J 2 r 


(4.2) 


In Fig. 4 we show how well the two methods above can extract fo: a result obtained from method¬ 
ologies that have different systematic errors is very robust. The same plot also displays a hyper¬ 
scaling fits of the data obtained from the Polyakov loop correlators. The fit has a very good y 2 /Ao\ 
and the fitted mass anomalous dimension y m = 0.3(1) is compatible with the rest of the spectrum. 
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Figure 4: Suing tension values in units of the lattice spacing for /3 = 4.0. Estimates obtained from Polyakov 
loop correlators and rectangular Wilson loop correlators are shown to be compatible for a range of volumes 
and masses. The Polyakov loops estimate is fitted with a hyperscaling form and gives y m = 0.3(1), broadly 
compatible with the rest of the hadronic spectrum. 


5. Conclusions 

Finite-size hyperscaling formulae describe the lattice data for the pseudoscalar mass and decay 
constant and for the vector mass, resulting in a universal value for the anomalous mass dimension. 
Both the flavor-singlet scalar mass and the string tension are compatible with these predictions. 
These are indications that the twelve-flavor SU(3) gauge theory simulated on the lattice corre¬ 
sponds to a conformal theory in the continuum and infinite volume limit. We also note that chiral 
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perturbation theory does not apply to this context for the reason that the expansion parameter is 
much larger than unity and that the would-be pions arc no longer the lightest degrees of freedom 
of the theory. 

Other theories that are thought to be close to or inside the conformal window seem to have 
a spectrum that contains a relatively light flavor-singlet scalar particle. If this feature continues 
to hold as the lattice studies become more precise and explore larger volumes at smaller fermion 
masses, it will be important to develop a fully fledged effective held theory that can take into 
account a scalar degree of freedom interacting together with the pseudoscalar ones. This is relevant 
for Walking Technicolor models, Composite Higgs models and any theory with light scalars. 
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